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An extended SU(3) shell model that for the first time explicitly includes unique-parity levels is 
introduced. Shell-model calculations for the isotopes of 64 Ge and 68 Se are performed where valence 
nucleons beyond the N=28=Z core occupy levels of the normal parity upper- fp shell (/5/2,P3/2, P1/2) 
and the unique parity g g / 2 intruder configuration. The levels of the upper- fp shell are handled within 
the framework of an m-scheme basis as well as its pseudo-SU(3) counterpart, and respectively, the 
g 9 / 2 as a single level and as a member for the complete gds shell. It is demonstrated that the 
extended SU(3) approach allows one to better probe the effects of deformation and to account for 
many key properties of the system by using a highly truncated model space. 



I. INTRODUCTION 

The nuclear shell model [l[ has been applied success- 
fully for the description of various aspects of nuclear 
structure, in large part because it is based on a minimum 
number of assumptions. Although direct diagonalization 
of the Hamiltonian matrix in the full Hilbert space would 
be desirable, the dimensionality of such a space is often 
too large to allow calculations of this type to be done. Re- 
cently, in order to relax this restriction dramatically, var- 
ious stochastic approaches, for instance, the Shell-Model 
Monte Carlo Method Q, have been suggested. Alterna- 
tively, algebraic models using the symmetry properties 
of the systems under investigation have been developed 
(e.g.||). 

Intruder levels are present in heavy deformed nuclei 
where the strong spin-orbit interaction destroys the un- 
derlying harmonic oscillator symmetry of the nuclear 
mean-field potential. The role they play for the over- 
all dynamics of the system has been the topic of many 
questions and debates 0, H, 0, @| ■ Until now, the prob- 
lem has been either approached within the framework 
of a truncation-free toy model ,5j or by just consider- 
ing the role of the single intruder level detached from 
its like-parity partners [H, 0] • It was argued in [f| that 
particles in these levels contribute in a complementary 
way to building the collectivity in nuclei. However, some 
mean-field theories suggest that these particles play the 
dominant role in inducing deformation 0]. In order to 
build a complete shell-model theory, these levels need to 
be included in the model space especially if experimen- 
tally observed high-spin or opposite-parity states are to 
be described. 

Until recently, SU(3) shell-model calculations - real 
SU(3) 3 for light nuclei and pscudo-SU(3) 0] for heavy 
nuclei - have been performed in either only one space 
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(protons and neutrons filling the same shell, e.g. the ds 
shell) or two spaces (protons and neutrons filling differ- 
ent shells, e.g. for rare-earth and actinide nuclei). Var- 
ious results for low-energy features, like energy spectra 
and electromagnetic transition strengths, have been pub- 
lished over the years [§, [l(| [H| ■ These applications con- 
firm that the SU(3) model works well for light nuclei and 
the pseudo-SU(3) scheme, under an appropriate set of 
assumptions, for rare-earth and actinide species. Up to 
now, SU(3)-based methodologies have not been applied 
to nuclei with mass numbers A = 56 to A = 100, which is 
an intermediate region where conventional wisdom sug- 
gests the break down of the assumptions that underpin 
their use in the other domains. In particular, the gg/ 2 
intruder level that penetrates down from the shell above 
due to the strong spin-orbit splitting appears to be as 
spectroscopically relevant to the overall dynamics as the 
normal-parity /5/2,P3/2iPi/2 levels. Specifically, in this 
region the effect of the intruder level cannot be ignored 
or mimicked through a "renormalization" of the normal- 
parity dynamics which is how it has been handled to date. 



The upper- fp + gg/ 2 shell is the lightest region of nu- 
clei where the intruder level must be taken into account. 
Its presence poses a significant challenge and opens a se- 
quence of questions, many of them still unanswered. For 
example, if the pscudo-SU(3) symmetry proves to be a 
good scheme for characterizing upper- fp shell configura- 
tions, should one integrate the gg / 2 intruder into this pic- 
ture by treating it as a single j-shell that is independent 
of couplings to the other members (37/2, ^5/2; ^3/21 S1/2) 
of the gds shell, or should one take the complete gds 
shell into account? While the treatment of the gg/ 2 in- 
truder as a single additional orbit is possible for only a 
few exceptional nuclei within the m-scheme shell-model 
calculations, an SU(3) approach to the problem could al- 
low the inclusion of the gg / 2 intruder as a member of the 
full gds shell and can be applied to all upper-fp nuclei. 
It is the purpose of this work to introduce and establish 
the benefits of a new and extended SU(3) shell model 
which, for the first time, explicitly includes particles from 
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FIG. 1: (Color online) Single-particle occupation numbers for eigenstates of the g.s. band of 64 Ge calculated in different 
restricted model spaces (from (a) to (e)) and the full space (f). The labels TPN over each restricted-space calculation represent 
the maximum total number of particles T, and the maximum number of protons P (and neutrons N) allowed in the intruder 
g g/2 level. 



the complete unique-parity sector and therefore can be 
used to explore the role that intruder levels play in the 
dynamics of the system. Calculations are performed for 
two nuclei which are of major importance in astrophysics, 
namely, the waiting-point N = Z nuclei 64 Ge and 68 Se 
fl2| . In addition, 68 Se is known to be among the nuclei 
for which shape coexistence effects have been reported 
O [3] ■ Both the strengths and the limitations of the 
model are demonstrated and discussed. 

First, a motivation for our approach is presented and 
the appropriate choice of model space is established by 
calculating, using a realistic Hamiltonian, the occupan- 
cies of the single-particle orbitals and the symmetry prop- 
erties of these typical uppcr-fp shell nuclei. Then, fol- 
lowing a brief introduction of the basics of the extended 
SU(3) model, calculations for the energy spectra, B(E2) 
transition strengths and the wave function content are 
performed and results are compared with the realistic 
predictions. Finally, discussion about various future ap- 
plications of the model is included. 

II. A REASONABLE APPROACH FOR THE 
DESCRIPTION OF UPPER FP-SHELL NUCLEI 

To benchmark the benefit of the SU(3) scheme in this 
region (pseudo-SU(3) for the upper- fp shell and normal 



SU(3) for the g g / 2 configurations extended to the full 
gds shell), we first generated results in a standard m- 
scheme representation for both nuclei, 64 Ge and 68 Se, 
with the 8 (4 protons + 4 neutrons) and 12 (6 protons 
+ 6 neutrons) valence nucleons, respectively, distributed 
across the Pi/ 2 ,p 3 /2, J '5/2, 39/2 model space with the f 7/2 
level considered to be fully occupied and part of a 56 Ni 
core. The Hamiltonian we used is a G-matrix with a phe- 
nomenologically adjusted monopole part [TR HH ] that in 
many cases describes the experimental energies reason- 
ably well. Specifically, this upper- fp + gg/ 2 shell inter- 
action was succesfully used in t he p ast to obtain quite 
good results for nuclei like 62 Ga [17j, an d 76 Ge and 82 Se 
fla ]. Later, it was applied for exploring the pseudo-SU(4) 
symmetry in the region from the beginning of the u ppe r 
fp-shell up to N = 30 and for describing beta decays [l8| ■ 

Calculations with different cuts of the full model space 
were done in order to estimate the occupancy of the 
single-particle levels and thus to evaluate the relative 
importance of various configurations for describing es- 
sential nuclear characteristics. In Fig. [TJ a comparison 
is made between the occupancies for 64 Ge as determined 
in a restricted basis, where at most two particles (pro- 
tons/neutrons) are allowed in the gg/ 2 level and the full- 
space resuls. The upper (yellow) bars show the contri- 
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FIG. 2: (Color online) Pseudo-SU(3) content of the low-lying 
states in (a) the g.s. band of 64 Ge, (b) and (c) - the K — 2 + 
band of 64 Ge, and (d) the g.s. band of 68 Se using the renor- 
malized counterpart of the G-matrix realistic interaction. 



bution to occupations from basis states with an occu- 
pied intruder level while the lower (blue) portion rep- 
resents those where the intruder level is empty. The 
calculated results suggest that the occupancy probabil- 
ity for the intruder level is approximately 0.3 particles 
for the low-lying states of 64 Ge. Calculations with no 
particles allowed in the intruder level or with just one 
identical-particle (or proton-neutron) pair (Fig. [T](a) and 
(b)) cannot describe either its occupancy or the gradual 
change in the occupancy of the single-particle levels in the 
ground-state (g.s.) band that is found in the full- model- 
space results. However, using a restricted space with at 
most two identical particles occupying the intruder level 
(in Fig. HJc)) is sufficient to describe both features as 



well as the low-energy spectrum and the B(E2) transi- 
tion strengths . Similar results were observed for the 
K = 2 + band of this nucleus. As expected, calculations 
for 68 Se performed in a truncated basis with at most 2 
nucleons allowed in the intruder level produce a slighlty 
higher value of the gg/ 2 occupancy compared to Ge. 

Next, the goodness of the pseudo-SU(3) symmetry in 
these nuclei was tested using a renormalized version of 
the same realistic interaction in the pf 5 / 2 space fl8j . The 
matrix of the second-order Casimir operator of pscudo- 
SU(3), C2 = j(3L 2 + Q-Q), was generated and the 
method of moments [2fj| used to diagonalize this matrix 
by starting the Lanzcos procedure with specific eigen- 
vectors of the Hamiltonian for which a pseudo-SU(3) de- 
composition was desired. Although the procedure pro- 
vides distributions only over the Gi values (and not over 
the actual (A, fj,) irreducible representations (irreps)) this 
analysis is quite useful and gives valuable information 
about the structure of certain eigenstates. 

The distribution of the second order Casimir operator 
C2 of pseudo-SU(3) yields contributions of about 50-60% 
from the leading pscudo-SU(3) irrcp in the g.s. band of 
64 Ge (Fig. HJa)) which suggests that the pseudo-SU(3) 
symmetry is quite good. In the K = 2 + band (Fig. [2] 
(b) and (c)) this contribution appears to be somewhat 
lower, ranging from approximately 37% for the state 
to about 62% in the 3^~ state. The analysis also reveals 
that using only five irreps which have the highest C2 value 
one may take account of at least 70% and up to about 
95% of the wavefunction for the states in these bands. 

In the case of 68 Se, the outcome turns out to be quite 
similar for the states from the g.s. band (Fig. [D^d)). 
Although the irreps with the maximal value of C2 = 180 
participate with only between about 40% and 50%, the 
first eleven irreps with distinct values of A and [i account 
for 88-93% of the wavefunction. In addition, the Oj 4 " state 
at 2.51 MeV is also dominated (64%) by irreps with the 
biggest C2 value. However, other states are predicted to 
be highly-mixed SU(3) configurations. This includes the 
O2 state found at 1.05 MeV - a value very similar to the 
ones reported in [3, EH for a low-lying state of prolate 
shape. A recent analysis reveales that many low-lying 
states in other N ~ Z nuclei also have good pseudo- 
SU(3) symmetry which further underscores the value of 
using symmetry-based truncation schemes [l9| . 

In summary, the outcome of the m-scheme study has 
demonsrtated that only part of the configurations are 
relevant for the structure of the low-lying states in the 
upper-fp shell nuclei. This along with the pseudo-SU(3) 
spectral decomposition of the states establishes the va- 
lidity of a SU(3)-bascd truncation scheme in N ~ Z 
upper-fp shell nuclei. As for the unique-parity spaces, 
the so-called quasi- SU(3) concept @, [Tfl, HH 

may also be 

applicable, however, this question is deferred to a follow- 
on investigation. Since (as we will see below) the role of 
the unique-parity spaces for the nuclei we deal with, is to 
introduce only some high-order effects, it is relatively safe 
to accept this to be true throughout the current work. 
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III. AN EXTENDED SU(3) MODEL WITH 
EXPLICITLY INCLUDED INTRUDER LEVELS 

Following the series of arguments and motivations pre- 
sented in the previous section, we can now introduce the 
basics of the extended SU(3) shell model. Like its early 
precursors 0, E|, it is also a microscopic theory in the 
sense that both SU(3) generators - the angular momen- 
tum (Lp, fi — 0, ±1) and quadrupole (Q^; [i — 0, ±1, ±2) 
operators - are given in terms of individual nucleon coor- 
dinate and momentum variables. However, the model 
space has a more complicated structure than the one 
used in earlier models based on the SU(3) symmetry. 
Specifically, it consists of two parts for each particle 
type, a normal (N) parity pseudo-shell {fs/2,P3/2,Pi/2 
— > d 5 /2, ^3/2 j S1/2) and a unique (or abnormal) (U) par- 
ity shell composed of all levels of opposite parity from 
the gds shell above. (Since the normal-unique space dis- 
tinction will be obvious from context, we will not place 
tildes over pseudo-space labels as is normally done.) 

The many-particle basis states 



where fku 



is essential and its contribu- 



tor; a ! ,}p(A, p)kL, {SV, S„}S; JM) 



(1) 



are built as SU(3) proton (it) and neutron (v) cou- 
pled configurations with well-defined particle num- 
ber and good total angular momentum. Here, the 
proton and neutron quantum numbers are indicated 
by a a = {a aN , a aU }p a {\ a , p a ), where the a aT = 
N aT [f aT }a aT (X (JT , p aT ) are the basis-state labels for the 
four spaces in the model (a stands for ir or is, and r 
stands for N or U). In the last expression, N aT denotes 
the number of particles in the corresponding space, [/ CTT ] 
- the spatial symmetry label and (A<rr, /•*<"-) - the SU(3) 
irrep label. Multiplicity indices a aT and p a count dif- 
ferent occurences of (Ao-tjAW) in [far] and in the prod- 
uct {(Ao-jvjMo-jv) x (XcrU,Vcru)} -* (Ao-./v), respectively. 
First, the particles from the normal and the unique spaces 
are coupled for both protons and neutrons. Then, the re- 
sulting proton and neutron irreps are coupled to a total 
final set of irreps. The total angular momentum J results 
from the coupling of the total orbital angular momentum 
L with the total spin S. The p and k are, respectively, the 
multiplicity indices for the different occurences of (A, p) 
in {(X n ,fi n ) x (X v ,fj, v )} and L in (A, p). 
The Hamiltonian 



H = V(iC" -GS* T iS ar 



<7. T 



2 



Q.Q : +a,J 2 + bKj 



- G(^2 S aT ^S aT ' +Y,S™' Tt S™ y ) ( 2 ) 

a.r^r' t,t' 

includes spherical Nilsson single-particle energies 

H sp =J2(Ho + -SW + ArrlL) ( 3 ) 

as well as the quadrupole-quadrupole and pairing inter- 
actions (within a shell and between shells) plus two rotor- 
like terms that arc diagonal in the SU(3) basis. In gen- 
eral, the harmonic oscillator term, Hq = Hui(r]i tTT + |) 



^41/3 22| j 

tion does not cancel out when more than one possi- 
ble distribution of particles over the shells is involved. 
The colons in the quadrupole operator notation repre- 
sent normal-ordered operator since all the one-body ef- 
fects have already been taken into account by the single- 
particle terms in the Hamiltonian. In addition, in first 
approximation, the quadrupole operator in the normal- 
parity spaces is related to its pseudo counterpart by 
Qo-n ~ ^jf-QaN with i) equal to 2 for both protons 
and neutrons and Q = Q nN + Q nU + Q„ N + Q uU f» 

1.5 QitN + QttIJ + 1-5 QvN + QvU- 



TABLE I: Parameters (in MeV) used in the extended SU(3) 
model Hamiltonian. 



Nucleus 


G 


X 


a 


b 


64 Ge 


0.280 


0.0176 


-0.002 


0.020 


68 Se 


0.263 


0.0152 


-0.002 


0.000 



The second line in Eq. @ consists of pairing terms that 
are included for the first time in SU(3) shell-model calcu- 
lations. In particular, the first term represents the scat- 
tering of an identical-particle pair between the normal- 
and unique-parity spaces. The second one stands for 
the proton-neutron pairing (or simply pn-pairing) inter- 
action within the normal- or unique-parity space (terms 
with r = t') and for the pn-pair scattering between the 
normal- and unique-parity spaces (terms with r 7^ r'). 
The formulae used in this work for each type of pairing 
operator can be found in the Appendix [SJ It is worth 
mentioning that, contrary to many calculations involv- 
ing the pairing interaction, these expressions are exact. 
Finally, the two rotor-like terms J 2 and Kj (the square 
of the total angular momentum and its projection on the 
intrinsic body-fixed axis) are used to "fine tune" the en- 
ergy spectra, adjusting the moment of inertia of the g.s. 
band and the position of the K = 2 + bandhead, respec- 
tively. Their strengths are the only two parameters fitted 
in this work. 

The single-particle terms together with the proton, 
neutron and proton-neutron pairing interactions mix the 
SU(3) basis states, which allows for a realistic description 
of the energy spectra of the nuclei. The values of the pa- 
rameters used in Hamiltonian can be found in Table 
HJ The single-particle energies in the Hamiltonian for 
the normal spaces are fixed with the numbers provided 
by the uppcr-fp shell single-particle energies and for the 
strengths in the unique-parity spaces the numbers from 
systematics are used [22|. The values for the parameters 
G and x m the Hamiltonian which are taken from [23[ 
are found to be in agreement with the ones [H, [24[ used 
in previous calculations for some ds-shcll and rare-earth 
nuclei. For simplicity, we take both identical-particle and 
proton-neutron pairing strengths to be equal. 
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TABLE II: The irreps in the coupled proton-neutron model space for 64 Ge and 68 Se that were used in the extended SU(3) 
shell-model calculations. The subscripts for each spin value denote the multiplicity, that is, the number of different ways the 
corresponding irrep can be constructed. 



[N n N, N-n-u; NyN, Nvu] total (A, ^multiplicity 



[4,0; 4,0] 


(8,4)0! 


64 Ge 

(9,2)0i,l a 


(10,0)0 2 ,li,2i 


(6,5)0i,l 2 


(7,3)0 6 ,1 8 ,22 


[4,0; 2, 2] 
([2, 2; 4,0]) 


(16, 2)0 2 


(17,0)12 


(14,3)06,16 


(15, l)0i2,li 5 ,2 4 


(12,4)0i8,li 8 ,2 4 


[3,1; 3,1] 


(16, 2)02,13,2! 


(17,0)02,13,2! 


(14,3)0 8 ,li 5 ,2 5 


(15,l)0l6,l27,2l 


(12,4)0 28 ,l 4 3,2l 5 


[2, 2; 2. 21 

L 5 J J J 


(24, 0)0i 


(22, 1)0 3 , 1 4 


(20, 2)0ifi, lis, 2 4 


(21, 0)0q, lis, 2 4 


( lo. O (0,19 . lftJ. . Zl fi 


[6 0-6 01 


(12 0l0i 

(l,10)0i,l 2 

(2,8)0 8 ,li2,2 4 


68 Se 

CO 12)0i 

(6,6)0 4 ,l 3 ,2i 


(9 3)09 lo 
(7,4)0 6 ,lii,2 3 


(3 9)09 I9 
(4, 7)0 6 ,lii,2 3 


cio DO1 19 

(8, 2)0g,lx2,24 


[6,0; 4, 2] 
([4, 2; 6,0]) 


(18, 2)0 2 

(17,1)012,116,24 


(19,0), la 
(13,6)0io,li 6 ,2 4 


(15,5)0 2 ,l2 


(16, 3)0 8 ,li2,2 2 


(12,8)02 


[5,1; 5,1] 


(18,2)0 2 ,l3,2i 

(17, 1)0X8, l 30 ,2l4 


(19, 0)0 2 , la, 2i 

(13,6)0 2 0,l33,2l5 


(15,5)04,1 6 ,2 2 


(16, 3)0l6,l27,2l3 


(12,8)0 2 ,l 3 ,2i 


[4, 2; 4, 2] 


(24,4)0x 

(24,l)0l9,l 3 l,2l3 


(25,2)0x,l2 
(20,6)0i6,lis,2 e 


(26,0)0 2 ,li,2i 


(22,5)0 3 ,l4 


(23, 3)Oi 4 ,l 2 o,2 6 



IV. RESULTS AND DISCUSSION 



Calculations within the framework of the extended 
SU(3) model were performed using irreps from 5 types 
of configurations - for example, [N^n, N„jj'i N u n, N u jj] 
= [4,0; 4,0], [4,0; 2, 2], [2, 2; 4,0], [3,1; 3,1] and [2, 2; 2, 2] 
for the 64 Ge case. (See Table llll where the list of config- 
urations for 68 Se is also given.) For each of these groups, 
irreps in the proton and neutron spaces with (pscudo-) 
spin S aT = 0, 1/2, 1 and 3/2 in both the normal- and 
the unique-parity spaces were generated. Then, from all 
the possible couplings between these we chose those with 
the largest value of the second order Casimir operator of 
SU(3) and spin S = 0, 1 and 2. Here, we present results 
obtained with five (seven) coupled proton-neutron irreps 
with distinct values of A and fi for each distribution of 
particles between the normal and unique spaces for 64 Ge 
( 68 Se). (This number was even pushed up to eleven for 
the [6,0;6,0] configuration in 68 Se). The complete set, 
listed in Table [TTl consists of 492 (580) coupled irreps in 
the case of 64 Ge ( 68 Se). This number is bigger by a fac- 
tor of about 20 than the one typically handled up to now 
within the framework of the SU(3) model. Some of the 
coupled irreps can be constructed in more than one way. 



For example, the irrep (A, /j,)S = (10, 0)0 can be obtained 
by coupling the (A,r, nJ)S v X {\ u , n v )S v = (4, 2)0 x (4, 2)0 
or (5,0)1 x (5,0)1 proton and neutron irreps. 

For both 64 Ge and 68 Se, proton-neutron configurations 
with no particles in the unique space are found, as ex- 
pected, to lie lowest and determine, by-and-large, the 
structure of the low-lying cigenstatcs. Only a small por- 
tion of all proton-neutron coupled irreps - 27 (112) in the 
case of 64 Ge ( 68 Se) - belong to these types of configura- 
tions, which we will refer to as the dominant ones. Since 
the only possible irrep in the unique-parity spaces for 
this case is {Ku,^u) = (0, 0) (and {\ vU ,Hvu) = (0, 0)), 
these configurations are the exact pseudo analog of the 
ones encountered in the ds-shell nuclei 24 Mg and 28 Si 
which have been studied earlier Q . Using only the prin- 
cipal part of the Hamiltonian namely, the part with 
both rotor term strengths equal to zero, we are able to 
provide a good description of the low-lying states. Specif- 
ically, all the energies from the g.s. bands (with the ex- 
ception of the 2+ state in 68 Se) differ by no more than 
15% from the experimental values [25|. In order to con- 
form with this result and prevent any further changes in 
the structure of the wave function, the range of values for 
the parameters a and b were severely restricted so that 
these terms only introduce small ( "fine tuning" ) changes 
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of a similar quality is provided by the extended SU(3) 
model. The existence of two prolate bands, as predicted 
by the calculations with the realistic interactions, is also 
observed, that is, a g.s. K = + and an excited K = 
2 + band, both dominated by the (8, 4) irrep. The first 
excited + (0^~) state, not reported yet experimentally, 
is found at 2.39 MeV which is higher than the prediction 
made by the realistic interactions. 

Consistent with the outcome for 64 Ge, in the case of 
68 Se we found a reasonable description for the energies of 
the states from the g.s. band (Fig. Even the use of a 
restricted space with at most 2 nuclcons allowed in the in- 
truder c/g / 2 level produces result which reflects some basic 
characterisics of the full-space spectrum reported in [l3T ] . 
For example, the first excited + state (0^) is also posi- 
tioned below the 2^ state. A new feature observed in our 
results is that the 0;|~ state at 2.51 MeV was found to be 
dominated by the shapes with C2 = 180 (see Section II). 
Within the framework of the extended SU(3) shell model, 
68 Se is predicted to be a mid-shell nucleus, a fact which 
may explain the existence of shape coexistence effects. 
Unlike the case of 28 Si 0], now the g.s. band is domi- 
nated by the irrep (12,0) which corresponds to a prolate 
shape. This result mainly follows due to the presence of 



FIG. 3: (Color online) Low-energy spectra of 64 Ge obtained 
with (from left to right) the realistic interaction in the full 
P.f 5 /2Q9/2 and p/5/2 model spaces compared with experiment 
[23 ] and the extended-SU(3)-model results. 



to the overall fit. 

Proton-neutron configurations with two and four par- 
ticles in the unique-parity space prevail at higher ener- 
gies. The former starts to dominate from about 3.5 MeV 
(5 MeV) in the case of 64 Ge ( 68 Se), usually for states 
of higher spin values, and the latter at even higher en- 
ergies. The amount of mixing found between configu- 
rations with different distribution of particles is due to 
the pair-scattering interactions between the normal- and 
unique-parity spaces. While the expected behavior of 
the unique-parity-space occupancy is observed - it goes 
up with the rise of both G and/or \ strengths - with the 
choice of parameters from Table Q] the absolute values ap- 
pear to be underestimated by at least a factor of 3 in 64 Ge 
(see the result shown in Fig. [I} and by even more for the 
case of 68 Se. This may indicate that the model space 
has to be further expanded to acommodate more mixing 
from the pair-scattering interaction terms, that the pn- 
interaction indeed should enter with different, possibly 
bigger strength than the identical-particle pairing, or the 
possible need to include other terms in the Hamiltonian 
like the quadrupole and isoscalar pairing interactions. 

Results for the excitation spectra of 64 Ge are presented 
in Fig. [3J The realistic G-matrix interaction gives a 
reasonable result for the low-lying states consistent with 
the one obtained in [l3|, (23|. Moreover, a description 
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FIG. 4: (Color online) Low-energy spectra of 68 Se obtained 
with (from left to right) the realistic interaction in the re- 
stricted P/5/2S9/2 (at most 2 particles allowed in the intruder 
gg/2 level) and full p/5/2 model spaces compared with exper- 
iment (25| and the extended-SU(3)-model results. 
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FIG. 5: (Color online) Role of the pn-pairing and the pair- 
scattering terms for the states in the g.s. band of 64 Ge and 
68 Se: (a) both pair-scattering and pn-pairing contributions 
excluded, (b) only pair-scattering contribution excluded and 
(c) total interaction. 

the orbit-orbit terms in the Hamiltonian and is in agree- 
ment with some earlier discussions [H, [2?) . Specifically, 
it favors the scenario in which the lower eigenstates in 
the g.s. band are prolate and throughout the band the 
shape changes to oblate [26| . Because of the nature of the 
leading representation, the model can not easily account 
for a K = 2+ band with the same shape characteristic, 
neither can it give a simple explanation for a low-lying 
K = + band, facts which are in support of the realistic 
prediction made in Section II for a highly-mixed nature 
of the O2 as well as many other low-lying states in this 
nucleus. With only the J 2 term used in the Hamiltonian 
for adjusting the energies, the 0^" state is predicted by 
the model at 1.55 MeV. 

The effect of adding the proton-neutron pairing and 
the pair-scattering terms in the Hamiltonian is illustrated 
in Fig. \5\ For our choice of model space ( Table HI)) , the 
results for the g.s. bands demonstrate comparable size 
effects from both interactions, especially for 64 Ge, with 
the only clear exception being the and 6^ states in 
68 Se. When a smaller number of irreps is included in the 
calculation, the pn-pairing interaction has a much bigger 
impact, an effect mainly visible for higher-spin states, 
while the role of the pair-scattering terms is strongly di- 
minished. 

Electromagnetic transition strengths are normally cal- 
culated with the E2 transition operator of the form 

T{E2) « y/E/^AVSfaVZ-tlQvN + e„Q nU 
Vv + 1 ~ 

+ &v— QvN + e-vQvu) (4) 

nu 

Instead, in this work we simply used the single domi- 
nant component in the pseudo-SU(3) expansion of the 
quadrupole operators in the normal-parity space. The 
effective charges e n and e„ were taken as e T = 1.5 and 
e v = 0.5 for the two versions of the realistic interac- 



TABLE III: B(E2) transition strengths for 64 Ge in units of 
e 2 /m 4 calculated using the G-matrix interaction in full p/5/2 
and p/5/239/2 model spaces, and the extended SU(3) model. 
Entries in parentheses show the result when only the normal 
spaces are used in the calculations. 



(J + 2)+ -» J+ 


Pf 5/2 


P/5/2S9/2 


Ext. SU(3) 


2+ 0+ 
4+ -> 2+ 

6 + — > 4 + 

8+ -> 6+ 


257.22 
332.54 
340.51 
303.31 


253.91 
342.51 
356.92 
320.14 


292.80 (280.10) 
346.26 (334.10) 
380.39 (370.56) 
273.84 (268.08) 


4+ -> 2+ 
8+ -> 6+ 


89.26 
164.23 
92.12 


93.13 
144.19 
84.38 


67.25 (65.73) 
207.18 (204.78) 
74.79 (79.39) 


(J+l)+ -» J+ 




371.15 
238.48 
159.44 


357.79 
240.40 
161.24 


505.27 (493.39) 
137.48 (135.48) 
10.26 (10.17) 


T+ -> 7+ 


2+ -> 0+ 

2+ -> 2+ 

3+ -> 2+ 

4+ -> 4+ 

4+ -> 2+ 
^7 "g.s. 


1.98 
251.68 
4.21 
72.10 
18.86 


1.42 
241.41 
3.40 
74.69 
19.31 


5.71 (5.74) 
183.16 (178.96) 

9.90 (9.93) 
47.11 (46.89) 
6.70 (6.75) 



TABLE IV: B(E2) transition strengths for the states in the 
g.s. band of 68 Se in units of e 2 /m 4 calculated using the G- 
matrix interaction in full p/5/2 model space and the extended 
SU(3) model. Entries in parentheses show the result when 



only the normal spaces 


are used in 


the calculations. 


(J + 2)+ -> J+ 


P/5/2 


Ext. SU(3) 


2+ -> 0+ 

"g.s. u g.s. 


322.71 


354.17 (346.37) 


4+ -> 2+ 
^g.s. "g.s. 


448.07 


486.65 (477.18) 


6+ -> 4+ 


441.58 


473.89 (467.09) 



tion and the extended-SU(3) calculations. The overall 
agreement between the results for both nuclei (Table 
HTT1 for 64 Gc and Table [TV] for the g.s. band in 68 Se) 
is good, although some recent experimental findings for 
the 2+ s — > 0+ s transition strength in 64 Ge [2!| seem to 
be underestimated by approximately a factor of 1.4. The 
correct behavior of the intcrband transitions is also nicely 
reproduced. More significant deviations are observed for 
the transitions between members of the K = 2 + band 
and the (J + 1) + — > J + transitions in 64 Ge. These 
could be attributed to the fact that some of the states 
from this band (e.g. 4^~ and 6^) are found to be highly 
mixed with S = 1 irreps and differ more significantly 
from the rest thus displaying a less regular structure pat- 
tern throughout the band, for example, to what has been 
observed in the same bands of some rare-earth nuclei. It 
seems that the orbit-orbit interaction is the part of the 
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TABLE V: Model-space dimensions for the G-matrix calculations in the full (restricted) pj "5/259/2 space for 64 Ge ( 68 Se) as well 
as for the complete p/5/2 spaces and for the extended SU(3) shell model with the irreps listed in Table ITT1 (The entry marked 
with a * is smaller by about a factor of two due to our taking advantage of time-reversal symmetry, which we did not invoke 
in the other cases as machine storage for them was not an issue.) 
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p/5/239/2 


1,831,531 


1,728,929 


1,454,930 


1,090,581 


724, 318 
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24, 246 


14, 760 


6,183 


1,638 


Ext. SU(3) 


322 


1,421 

68 Se 


2,098 


2,225 


2,208 


P/5/239/2 


1,929,014* 


3,611,680 


2, 973, 404 


2, 138,391 




P/5/2 


93, 710 


81, 122 


52, 175 
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FIG. 6: (Color online) Wave- function decomposition of the 
calculated extended SU(3) eigenstates in the g.s. band of (a) 
64 Ge and (b) 68 Se. The leading irreps from the dominant type 
configurations are listed explicitly while the effect of those 
with less than 3% contribution for any state as well as from 
configurations with two and four particles in the unique space 
is represented with a dashed line. 



Hamiltonian responsible for this feature. Results when 
only the normal-parity spaces are included in the calcu- 
lation (shown in parentheses) reveal a contribution of the 
unique-parity sector of only up to 2-3%. An increase in 
this number is expected for higher-lying states or heav- 
ier nuclei where the dominant configurations are the ones 
with an occupied unique-parity space. 

Finally, let us look at the content of the eigenfunctions 
for the states in the different bands. In the g.s. (K = + ) 
band of 64 Ge, one can clearly see the dominance of the 
leading and most deformed SU (3) irrcp (8, 4) (Fig. [B^a)) 
which gradually declines throughout the band from about 
80% for J = 0+ s to less than 40% for J = 8+. Since the 
spin-orbit interaction is not as strong as in the case of the 
ds-shell nuclei, the mixing of irreps is smaller compared 
to the corresponding normal-SU(3) results for Mg and 
28 Si [{|. On the other hand, as mentioned above, the 
K = 2 + band follows a less regular pattern with some of 
the states being of highly mixed nature. The 0^ state is 
found to be dominated (85%) by the (A, /i), S = (9, 2), 1 
irrep. In the case of 68 Se, the leading irrep (12,0) con- 
tributes from 75 to 85% (Fig. [6jb)). A slight change 
in the type of Hamiltonian used may help establish the 
transition from states of prolate shape dominated by the 
irrep (12,0) in the g.s. band to ones where the (0,12) ir- 
rep prevails. To achieve this effect we need to add a term 
proportional to the third order Casimir invariant C3 of 
SU(3). It was demostrated earlier Q that this term is 
capable of adjusting the prolate-oblate band crossing by 
driving irreps with fi » A lower in energy than those 
with A >> (j,. The same term can also be used to fix the 
position of the first excited + state not assigned yet ex- 
perimentally but predicted by our G-matrix calculations 
to lie at 1.05 MeV. 

Although the extended-SU(3) calculations are per- 
formed in a model space that involves the whole gds shell, 
the basis is still much smaller in size even compared with 
the one used for realistic calculations in the pj 5/2*79/2 
space. This drastic reduction translates into the use of 
only hundreds or at most a few thousand basis states 
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(Tabic [V}. For example, the size of the basis used in the 
extended-SU(3) calculations for 64 Ge represents only be- 
tween 0.02% to 0.3% of that for unrestricted calculations 
in the pf 5/299/2 model space. This means that a space 
spanned by a set of extended-SU(3) basis states may be 
computationally manageable beyond the limit accessible 
for the modern full-space shell-model calculations as is 
the case for the combination of the uppcr-fp and the gds 
shells. 

While some refinements in the model certainly could be 
done (like trying different and more sophisticated types 
of Hamiltonians, using different strengths for identical- 
particle and pn-pairing interactions, etc.) and the role 
of the model-space truncation may be further explored, 
the results presented in this paper suffice to demonstrate 
that the SU(3) scheme in its extended formulation can be 
a valuable tool for studying nuclei of the upper-fp-shell 
region. Its benefits will show even more prominently in 
the more general and complicated case, namely, when 
the dominant configuration is no longer the one with an 
empty unique-parity space and in situations when two or 
more competing configurations are closer to one another 
in energy and as a consequence experience strong mixing. 

V. CONCLUSIONS 

In this paper we extended the usual pseudo-SU(3) shell 
model for upper-fp shell nuclei in two ways: firstly by 
integrating the gg/2 level into the dynamics; and sec- 
ondly by including the entire gds-shell organized via its 
SU(3) structure, which we dubbed the extended SU(3) 
shell model. While this work only deals with the sim- 
plest case in which one configuration (the one with no 
particles in the unique-parity space) dominates all others, 
it is still possible to appreciate the strength of this new 
approach. Specifically, the model offers a richer model 
space compared to the previous SU(3) schemes by tak- 
ing particles from the unique-parity space explicitly into 
account. As a result, the current approach presents an 
opportunity for a better description of the collectivity 
properties of the systems considered by reducing the ef- 
fective charge needed in the description of their B(E2) 
transition strengths. These results will be even more pro- 
nounced for heavier systems where the intruder space is 
expected to have higher occupancy. This approach also 
offers an opportunity to explore the role of the intruder 
levels in the dynamics of the system as in the current 
study they are treated on the same footing as the normal- 



parity orbitals. It is important to underscore that these 
advantages are accomplished within a highly truncated 
and symmetry-adapted basis, which possibly allows one 
to reach into otherwise computationally challenging (if 
not inaccessible) domains. 

The results for the nuclei 64 Ge and 68 Se, presented 
in this paper, demonstrate a close reproduction of vari- 
ous results obtained with a realistic interation. Specifi- 
cally, many of the states in the energy spectra and the 
B(E2) transition strengths are nicely reproduced. While 
the results are satisfactory for the states from the g.s. 
bands, there still seem to be some need for a more pre- 
cise description of the nuclear characteristics related to 
the properties of the eigenfunctions. These could be ad- 
dressed in the future by including some corrections with 
the use of more elaborated interactions. Nevertheless, the 
results certainly suggest that the extended SU(3) model 
can be a valuable tool in studying properties of nuclei of 
special interest from this region, such as those lying close 
to the proton drip line or/and actively participating in 
the processes of nucleosynthesys. They also point to an 
excellent opportunity to reveal the role the intruder lev- 
els play in the dynamics of the system in an exciting and 
completely new way, namely, considering their connec- 
tion to their like-parity partners within the framework of 
a severely-truncated symmetry-adapted model space. 

An extension of the SU(3) shell model for the rare- 
earth and actinide nuclei is also underway. It will be 
able to provide some valuable new information and a 
better understanding of the fragmentation and cluster- 
ization phenomena in the B(M1) transition strengths. 
Also, it will signifficantly reduce the values of the ef- 
fective charges used in estimates of the B(E2) transition 
strengths. In addition, an expected new emerging struc- 
ture of the states in the excited K = + bands could 
give an explanation for the enhanced B(E2) transition 
strengths to members of the g.s. band. Finally, the new 
model provides a powerful means of explanation for the 
abundance of low-lying K = + states found experimen- 
tally. 
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APPENDIX A: MONOPOLE PAIRING AND PAIR-SCATTERING IN THE FRAMEWORK OF THE 

EXTENDED SU(3) SHELL MODEL 



The expression for the identical-particle pairing within a shell (and the identical-particle pair-scattering between 
two shells) in terms of SU(3) irreducible operators is given in [30| : 

S^S"' = \ Yl EH'"'' VC21 + l)(2l' + 0)1; (r,, 0)i||(Ai, M i)10> ((0, rf)l'; (0, V ')l'\\ fa, A 2 )10) 

(Ai,/i 1 )(A 2 ,M2) " 
p(X,fi) 



x((Ai, M i)10;(M2,A 2 )10||(A,/i)10) p 



a L,o)i x a Lo)i 



(Ai,/ii)Si=0 



a (0,J)')| X a (0,r;')5 



(M2,A 2 )S 2 =0 



(Al) 

K= 1L=0 Mj=0 



where cr = ir or 1/ with ?y = r)' (r) ^ rf) and r = r' (r 7^ r') for the case of pairing (pair-scattering). Here, 
( ; || ) denotes a reduced SU(3) D SO (3) Clebsch-Gordan coefficient and a is a proper SU(3) tensor defined by 

p; ( \ri-\-j+Tn „ 

fl (0,i))(]m — a ()),0)/j-m- 

Using the fact that the 9 — A/i (9j) coefficients connect composite tensors corresponding to different coupling schemes 
of four SU(3) (SU(2)) tensors, one can derive the corresponding expressions for the pn-pairing (pair-scattering) 
operators. The final result is given by 



1 

4 



J2 £ (-) i - r V2£TT^V25TT(( f7 ,0)/;(0,7/)r||(A T ,^) K7r £) 

(A 7r ,^ 7r )(A l/ ,^ 1/ ) 
p'(A, M ) 

x ((77, 0)/; (0, t/^'II (A„, [i v )k v L) ((A w , y u 7r )K 7r £; (A^, ^)k^£|| (A, ^)10) p ' 



a (t? , )i X & (0,r,')| 



%,0)I x 6 (o,„')! 



p(A,/i),S=0; ,7=0 



K= 1L=0 M. T =0 



(A2) 



where a^(a) and b^{b) denote proton and neutron creation (annihilation) operators. Again, 77 = 7/ (77 7^ 77') and r = r' 
(t 7^ t') for pn-pairing (pn-pair scattering). The labels (A^,^) and (A„,/^„) represent the proton and the neutron 
SU(3)-coupled irreps and the symbols C and S stand for the orbital and spin angular momentum in the coupled 
proton (and neutron) spaces. 



[1] M. G. Meyer, Phys. Rev. 75, 1969 (1949); 78, 16 (1950); 

78, 22 (1950); O. Haxel, J. H. D. Jensen, and H. E. Suess, 

ibid. 75, 1766 (1949). 
[2] S. E. Koonin, D. J. Dean, and K. Langanke, Phys. Repts. 

577, 1 (1996). 

[3] J. P. Elliott, Proc. Roy. Soc. London, Ser. A 245, 128 

(1958); A 245, 562 (1958). 
[4] R. D. Ratna Raju, J. P. Draayer, and K. T. Hecht, Nucl. 

Phys. A202, 433 (1973). 
[5] J. Escher, J. P. Draayer, and A. Faessler, Nucl. Phys. 

A586, 73 (1995). 
[6] K. H. Bhatt, C. W. Nestor, Jr., and S. Raman, Phys. 

Rev. C 46, 164 (1992). 
[7] K. H. Bhatt, S. Kahane, and S. Raman, Phys. Rev. C 

61, 034317 (2000). 
[8] S. Aberg, H. Flocard, and W. Nazarewicz, Annu. Rev. 

Nucl. Part. Sci. 40, 469 (1990). 
[9] C. Vargas, J. G. Hirsch, and J. P. Draayer, Nucl. Phys. 

A690, 409 (2001). 
[10] C. E. Vargas, J. G. Hirsch, and J. P. Draayer, Nucl. Phys. 



A697, 655 (2002). 
[11] G. Popa, J. G. Hirsch, and J. P. Draayer, Phys. Rev. C 

62, 064313 (2000); C. Vargas, J. G. Hirsch, T. Beuschel, 

and J. P. Draayer, Phys. Rev. C 61, 031301(R) (2000). 
[12] H. Schatz et al., Phys. Rep. 294, 167 (1998); J. A. Clark 

et al., Phys. Rev. C 75, 032801(R) (2007). 
[13] K. Kaneko, M. Hasegawa, and T. Mizusaki, Phys. Rev. 

C 70, 051301(R) (2004). 
[14] Y. Sun, M. Wiescher, A. Aprahamian, J. Fisker, Nucl. 

Phys. A758, 765c (2005). 
[15] E. Caurier, F. Nowacki, A. Poves, and J. Retamosa, 

Phys. Rev. Lett. 77, 1954 (1996). 
[16] A. Abzouzi, E. Caurier, and A. P. Zuker, Phys. Rev. 

Lett. 66, 1134 (1991). F. Nowacki, Ph. D. thesis, ULP 

Strasbourg, 1995 (unpublished). 
[17] S. M. Vincent et al., Phys. Lett. B 437, 264 (1998). 
[18] P. Van Isacker, O. Juillet, and F. Nowacki, Phys. Rev. 

Lett. 82, 2060 (1999). 
[19] K. P. Drumev et al., to be submitted to Nucl. Phys. A. 
[20] R.R. Whitehead and A. Watt, J. Phys. G 4, 835 (1978); 



11 



R.R. Whitehead et al., in Theory and Application of mo- 
ment methods in Many-Fermion systems, edited by B. J. 
Dalton et al. (Plenum, New York, 1980); E. Caurier, A. 
Poves, and A. P. Zuker, Phys. Lett. B252, 13 (1990); 
Phys. Rev. Lett. 74, 1517 (1995). 

[21] A. P. Zuker, J. Retamosa, A. Poves, and E. Caurier, 
Phys. Rev. C 52, R1741 (1995); G. Martinez-Pinedo, A. 
P. Zuker, A. Poves, and E. Caurier, Phys. Rev. C 55, 
187 (1997); C. Vargas, J. G. Hirsch, P. O. Hess, and J. 
P. Draayer, Phys. Rev. C 58, 1488 (1998). 

[22] P. Ring and P. Schuck, The Nuclear Many-Body Problem 
(Springer, Berlin 1979). 

[23] K. Kaneko, M. Hasegawa, and T. Mizusaki, Phys. Rev. 
C 66, 051306(R) (2002). 

[24] M. Dufour, A. P. Zuker, Phys. Rev. C 54, 1641 (1996). 



[25] National Nuclear Data Center, 

|http://www.nndc. bnl.gov 
[26] K. C. Tripathy and R. Sahu, Int. Journ. Mod. Phys. E 

11, 531 (2002). 

[27] S. M. Fischer, D. P. Balamuth, P. A. Hausladen, C. J. 

Lister, M. P. Carpenter, D. Seweryniak, and J. Schwartz, 

Phys. Rev. Lett. 84, 4064 (2000). 
[28] J. P. Draayer and K. J. Weeks, Ann. Phys. (N.Y.) 156, 

41 (1984); O. Castanos, J. P. Draayer, and Y. Leschber, 

ibid. 180, 290 (1987). 
[29] K. Starosta et al., Phys. Rev. Lett. 99, 042503 (2007). 
[30] C. Bahri, J. Escher, and J. P. Draayer, Nucl. Phys. A592, 

171 (1995). 



